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Abstract 

In this paper, we present a framework for automatic generation of CHR solvers given the 
logical specification of the constraints. This approach takes advantage of the power of 
tabled resolution for constraint logic programming, in order to check the validity of the 
rules. Compared to previous wor k l|Apt and Monfroy 1999||Ringeissen and Monfroy 2000| 
[Abdennadher and Rigotti 2000| |XM€nnaHTi€r^m^^Ugotti20^^| l , where different meth- 
ods for automatic generation of constraint solvers have been proposed, our approach en- 
ables the generation of more expressive rules (even recursive and splitting rules) that can 
be used directly as CHR solvers. 

KEYWORDS: Rule-based constraint solver, automatic solver generation, tabled resolu- 
tion, constraint logic programming 



1 Introduction 

Constraint Handling Rules (CHR) l|Friihwirth 1998|l is a high-level language espe- 
cially designed for writing constraint solvers. CHR is essentially a committed-choice 
language consisting of multi-headed rules that transform constraints into simpler 
ones until they are solved. CHR defines both simplification of and propagation over 
user-defined constraints. Simplification replaces constraints by simpler constraints 
while preserving logical equivalence. Propagation adds new constraints, which are 
logically redundant but may cause further simplifications. Consider the constraint 
mm, where minlX, Y, Z) means that Z is the minimum of X and Y. Then typical 
CHR rules for this constraint are: 

min{X,Y,Z), Y<X ^ Z=Y, Y<X. 
min{X,Y,Z), X<Y ^ Z=X, X<Y. 
min{X,Y,Z) ^ Z<X, Z<Y. 
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The first two rules are simplification rules, while the third one is a propagation 
rule. The first two rules correspond to the usual definition of min. The first rule 
states that min{X,Y,Z), Y<X can be replaced by Z~Y, Y<X. The second has 
an analogous reading, and the third rule states that if we have mm(X, y, Z) then 
we can add Z<X, Z<Y to the current constraints. 

If such rules are in general easy to read, in many cases it remains a hard task to 
find them when one wants to write a constraint solver. Thus, several methods have 
been proposed to automatically generate rule-based solvers for constraints given 
their logical specification ( |Apt and Monfroy 1999, .Ringeissen and Monfroy 2000| 
[Abdennadher and Rigotti 2000||Xbdennadher and Rigotti 200 la] ). These approaches 

can help to find more easily interesting rules, and it has also been shown in ( [Abdennadher and Rigotti 2001b| ) 
that the rules produced automatically can lead to more efficient constraint reason- 
ing than rules found by programmers. 

In this paper, we propose a new method to generate CHR propagation and sim- 
plification rules. This work extends the previous rule-based solver generation tech- 
niques described in I jApt and~Monfroy 1999,.Ringeissen and Monfroy 2000l [Abdennadher and Rigotti 2000| 
[Abdennadher and Rigotti 2001 at . It allows to obtain more general forms of rules, 
even for constraints defined intensionally over infinite domains. 

The intuitive principle of the generation is the following. Consider that a solver 
S for some constraints, called primitive constraints, is already available. Other con- 
straint predicates, called user-defined constraints, are given and their semantics is 
specified by mean of a constraint logic program P (i.e., the constraint predicates 
are defined by clauses of P). Then we want to obtain CHR propagation and sim- 
plification rules for the user-defined constraints to extend the existing solver. The 
basic idea of our approach relies on the following observation: a rule of the form 
C ^ D is valid if the execution of the goal C, -^{D) finitely fails with program P 
and solver S. For the execution of such goals, we will use a tabled resolution for 
constraint logic programming IjCui and Warren 2000|l that terminates more often 
than execution based on SLD-like resolutions. 

We present three algorithms that can be integrated to build an environment to 
help developers to write CHR rule-based constraint solvers. Two of the algorithms 
focus on how to generate propagation rules for constraints given their logical spec- 
ification. The first algorithm generates only primitive propagation rules (i.e., rules 
with right hand side consisting of primitive constraints). The second algorithm 
extends the first one to generate more general propagation rules with right hand 
side consisting of both primitive and user-defined constraints. We also show that a 
slight modification of this algorithm allows to generate the so-called splitting rules 
(rule having a disjunction in their right hand side) supported by the extension of 
CHR called CHR^ ([Abdennadher and Schiitz 1998|l . The third algorithm focuses 
on transforming propagation rules into simplification rules to improve the time and 
space behavior of constraint solving. 

Organization of the paper. In Section 13 we present an algorithm to generate prim- 
itive propagation rules. In Section O we describe how to modify the algorithm 
to generate more general propagation rules. Section 0] presents a transformation 
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method of propagation rules into simplification rules. For clarity reasons, we will 
use an abstract representation for the generated rules. In Section |5l we present 
how these rules can be encoded in CHR. We discuss related work in Sectional and 
finally we conclude with a summary and possibilities of further improvements. 

2 Generation of Primitive Propagation Rules 

We assume some familiarity with constraint logic programming (CLP) IjJaffar and Maher 19941 
[Marriott and Stuckey 1998| ). We consider two classes of constraints, primitive con- 
straints and user-defined constraints. Primitive constraints are those constraints 
defined by a constraint theory CT and for which solvers are already available. In 
the following, we do not expect that the solver for primitive constraints is complete. 
User-defined constraints are those constraints defined by a constraint logic program 
P and for which we want to generate solvers. We assume that the set of primitive 
constraints is closed under negation, in the sense that the negation of each primitive 
constraint must be also a primitive constraint, e.g. = and ^ or < and >. In the 
following, we denote the negation of a primitive constraint c by not{c). 
In the rest of this paper, we use the following terminology. 

Definition 2.1 

A constraint logic program is a set of clauses of the form 

h < 5i, . . . , bji^ Ci , . . . , Cjri 

where /i, 6i, . . . , 6„ are user-defined constraints and ci, . . . , Cm are primitive con- 
straints, h is called left hand side of the clause. A goal is a set of primitive and 
user-defined constraints. An answer is a set of primitive constraints. The logical 
semantics of a constraint logic program P is its Clark's completion and is denoted 
by P* . A user-defined constraint is defined in a constraint logic program if it occurs 
in the left hand side of a clause. 

Definition 2.2 

A primitive propagation rule is a rule of the form Ci ^ C2 or of the form Ci => 
false, where C'l is a set of primitive and user-defined constraints, while C'2 consists 
only of primitive constraints. Ci is called the left hand side of the rule (Ihs) and C2 
its right hand side (rhs). A rule of the form Ci false is called failure rule. 

In the following we use the notation El((/)) to denote the existential closure of (p 
and 3_v(</') to denote the existential closure of (j> except for the variables in the set 
V. 

Definition 2.3 

A primitive propagation rule {di, . . . , (i„} {ci, . . . , Cm} is valid with respect to 
the constraint theory CT and the program P if and only if P* , CT \= IX^di — > 
3_v(Aj Cj), where V is the set of variables appearing in {c?i, . . . , A failure rule 
{di, . . . , (i„} => false is valid with respect to CT and P if and only if P*, CT |= 

-3(A,rf.). 

We now give an algorithm to generate valid primitive propagation rules. 
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2.1 The PRIM-MINER Algorithm 

The PRIM-MINER algorithm takes as input the program P defining the user- 
defined constraints. To specify the syntactic form of the rules, the algorithm needs 
also as input two sets of primitive and user-defined constraints denoted by Baseihs 
and Candihs , and a set containing only primitive constraints denoted by Candrts ■ 
The constraints occurring in Baseihs are the common part that must appear in 
the Ihs of all rules, Candihs indicates candidate constraints used in conjunction 
with Baseihs to form the Ihs, and Candrhs are the candidate constraints that may 
appear in the rhs. 

Note that a syntactic analysis of the constraint logic program P can suggest 
functors and constraint predicates to be used to form candidate constraints. 

The algorithm prim-miner is presented in Figure ITTl and generates a set of valid 
rules of the form C d where d is a single primitive constraint. Note that rules 
with the same Ihs, e.g. C => di, C ^ d2, ■ ■ ■ , C ^ dn, can be replaced by the single 
rule C => di,d2, . . . ,dn. 

The basic idea of the algorithm relies on the following observation: to be able 
to generate a failure rule of the form C false, we can simply check that the 
execution of the goal C finitely fails. Furthermore, while these rules are useful to 
detect inconsistencies, it is in general more interesting to propagate earlier some 
information that can be used for constraint solving, instead of waiting until a con- 
junction of constraints becomes inconsistent. Thus, for each possible Ihs C (i.e., 
each subset of Baseihs U Candihs) the algorithm distinguishes two cases: 

1. PRIM-MINER uses a CLP system to evaluate the goal C. If the goal finitely 
fails, then the failure rule C false is generated. 

2. Otherwise the negation of each candidate constraint d from Candrhs is added 
in turn to C and the goal C U {not{d)}) is evaluated. If the goal finitely fails, 
then the rule C ^ {d} is generated. 

In practice these goal evaluations are made using a bounded depth resolution 
procedure to avoid non-termination of the whole generation algorithm. 

The algorithm prim-miner uses a basic ordering to prune the search space and to 
avoid the generation of many uninteresting rules. This pruning relies simply on the 
following observation. If Ci ^ false is valid, then rules of the form C2 false, 
where Ci C C2 are also valid but useless. So the algorithm considers first the 
smallest Ihs with respect to set inclusion, and when it finds a valid failure rule 
Ci => false it discards from the Ihs candidates any C2 that is superset of Ci. 

At first glance, the procedure used to evaluate the goals issued by the algorithm 
may be considered as a classical depth-first, left-to-right CLP resolution. However, 
we will show in Section 12.21 and Section |21 that a tabled CLP resolution extends 
greatly the class of rules that can be generated, by allowing termination of the 
evaluation in many interesting cases. Additionally, it should be noticed that the 
execution on the underlying CLP system is not required to enumerate all answers 
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begin 

TZ the resulting rule set is initialized to the empty set. 
L is a list of all subsets of Candihs, 
in an order compatible with the subset partial ordering 
(i.e., for all Ci in L if C2 is after Ci in L then C2 (Z^ Ci). 

while L is not empty do 

Remove from L its first element denoted Cihs- 
if the goal {Baseihs U Cihs) fails 

with respect to the constraint logic program P then 
add the failure rule {Basems U Cihs => false) to TZ 
and remove from L all supersets of Cihs- 

else 

Let rhs be initialized to the empty set. 
for all d G Candrhs 

if the goal {Baseihs U Cihs U {not{d)}) fails 

with respect to the constraint logic program P then 
add d to the set rhs. 
endif 
endfor 

if rhs is not empty then add the rule {Baseihs U Cihs => rhs) to TZ endif 
endif 
endwhile 

output TZ. 

end 



Fig. 1. The prim- miner Algorithm 

since PRIM-MINER only performs a fail/succeed test, and thus the CLP system can 
stop after a first answer has been found. 

Example 2.1 

Consider the following constraint logic program which implements the predicate 
min. min{X, Y, Z) means that Z is the minimum of X and Y: 

min{X,Y,Z) ^ X<Y, Z=X. 
min{X,Y,Z) ^ Y<X, Z=Y. 

where < and = are primitive constraints with the usual meaning as total order and 
syntactic equality. 

The algorithm with the appropriate input generates (among others) the rule 

min{X,Y,Z), Y<X => Z=Y. 

after having checked that the execution of the goal min{X, Y, Z), Y<X, Z^Y fails 
by constructing the following derivation tree: 
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min{X,Y,Z), Y<X, Z^Y 





X<Y, Z=X, Y<X, Z^Y 



Y<X, Z=Y, Y<X, Z^Y 



false 



false 



Note that we assume that the constraint solver for < and = is able to detect such 
inconsistencies. 

Soundness and Completeness. The prim-miner algorithm attempts to extract all 
valid primitive propagation rules of the form Ci => {d} or Ci =^ false such that 
Baseihs C Ci, Ci \ Baseihs C Candihs, d G Candrhs, and there is no other more 
general failure rule (i.e., no valid rule C2 ^ false where C2 C Ci). In general, the 
algorithm cannot be complete, since the evaluation of some goals corresponding to 
valid rules may be non-terminating. In fact, this completeness can be achieved if 
more restricted classes of constraint logic programs are used to give the semantics 
of user-defined constraints and if the solver for the primitive constraints used by 
the underlying CLP system is satisfaction complete. 

The soundness of the algorithm (i.e., only valid rules are generated) is guaranteed 
by the nice properties of standard CLP schemes (|,Taffar and Maher 1994|l and tabled 
CLP schemes l|Cui and Warren 2000|l . An important practical aspect, is that even 
a partial resolution procedure (e.g., bounded depth evaluation) or the use of an 
incomplete solver by the CLP system, does not compromise the validity of the 
rules generated. 

2.2 Advantage of Tabled Resolution for Rule Generation 

Termination of the evaluation of (constraint) logic programs has received a lot of 
attention. A very powerful and elegant approach based on tabled resolution has been 
developed, first for logic programming (e.g., IjTamaki and Sato 1986IIWarren 1992|l ') 
and further extended in the context of CLP (e.g., IjCui and Warren 2000|l ). 

The intuitive basic principle of tabled resolution is the following. Each new sub- 
goal S is compared to the previous intermediate subgoals (not necessarily in the 
same branch of the resolution tree). If there is a previous subgoal / which is equiv- 
alent to S or more general than 5, then no more unfolding is performed on S and 
answers for S are selected among the answers of /. This process is repeated for all 
subsequent computed answers that correspond to the subgoal /. 

The use of such technique has not been widely accepted since if this leads to ter- 
mination in many more cases than execution based on SLD-resolution, this should 
be paid by some execution overhead in general. When using the algorithm PRIM- 
MINER we can accept a slight decrease of performance (since the solver is constructed 
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once) if this gives rise to an improvement of the termination capability which en- 
ables the generation of additional rules. Thus, tabled CLP can find very interesting 
applications in constraint solver synthesis. This will be illustrated in the following 
example. 

Example 2.2 

Consider the well-known ternary append predicate for lists, which holds if its third 
argument is a concatenation of the first and the second argument. 

appendix, Y,Z) ^ X=[], Y=Z. 

appendix, Y,Z) ^ X=[H\X1], Z=[H\Z1], appendix 1,Y, Zl). 
Let the input parameters of prim-miner be 
Baseihs = {appendix, Y,Z)} 

Candihs = {X^[],Y^[],Z=[],X^Y,X^Z,Y^Z,X^Y,X^Z,Y^Z] 

Candrhs = Candihs- 

Then the algorithm generates (among others) the following primitive propagation 
rule: 

appendix, Y,Z), Y=W X=Z 

by executing the goal appendix, Y, Z), Y~W, X^Z with a tabled CLP resolution. 
For a classical CLP scheme the resolution tree will be infinite, while in case of a 
tabled resolution it can be sketched as follows: 

appendix, Y, Z), Y^W,X^Z 




X ^ \\,Y^Z,Y^W,X^Z X^[H\Xl], Z=[H\Zl],appendiXl,Y, Zl),Y^W, X^Z 

false false 

The initial goal Gi = iappendiX ,Y , Z) , Y=W, X^Z) is more general than the 
subgoal G2 = iX={U\X\\, Z=[H\Z1], appendix 1, Y, Zl), Y=[], X^Z), in the 
sense that appendiA,B,C), D^[E\A], F=[E\C], B=[], D^F entails appendiA, B, 
C), B—W, A^C. So no unfolding is made on G2, and the process waits for answers of 
Gi to compute answers of G2. Since Gi has no further possibility of having answers, 
then G2 fails and thus Gi also fails. We refer the reader to IjCui and Warren 2000|l 
for a more detailed presentation of goal evaluation using a tabled CLP resolution. 

Using this kind of resolution prim-miner is also able to produce rules such as: 

appendix, Y,Z), X=Z Y=[]. 
appendix, Y,Z), Y^W X^Z. 
appendix, Y,Z), X^[] => Z^[]. 



We have run all examples presented in this paper, using our own implementation 
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of tabled CLP resolution according to the description of IjCui and Warren 2000|l . 
The running prototype is implemented in SICStus Prolog 3.7.1. 

2. 3 Selection of Interesting Rules and Performance of Generation 

Interesting rules. If we use the algorithm prim-miner as presented in Figure ITTl 
we obtain a set of rules that is highly redundant and thus not suitable neither 
for human-reading nor for its direct use as an executable constraint solver. This 
problem encountered in propagation rule generation has already been pointed out 
in ( [Abdennadher and Rigotti 2000| ). It is taken into account by applying the fol- 
lowing simplification technique on the set of rules produced by prim-miner: 

• The rules generated by PRIM-MINER are ordered in a list L using any total or- 
dering on the rule Ihs compatible with the 0-subsumption ordering IjPlotkin 1970|l 
(i.e., a rule having a more general Ihs is placed before a rule with a more spe- 
cialized Ihs). 

• Let 5 be a set of rules initialized to the empty set. For each rule Ci C2 in 
I (taken according to the list ordering) the constraint C2 is simplified to an 
equivalent constraint Csimp by the already known solver for Candrhs and by 
the rules in S. If Csimp is empty then the rule can be discarded, else add the 
rule Ci Csimp to S. 

• Output the set S containing the simplified set of rules. 

For example, using this process the set of rules 

{p{X) ^ r[X), p{X), q(X) ^ r(X), s(X, Y) => X=Y, X=a, Y=a} 
is simplified to 

{p{X) r{X), s{X, Y) ^ X^a, Y^a}. 

For clarity reasons this simplification step is presented separately from the al- 
gorithm PRIM-MINER, but it can be incorporated in prim-miner and performed 
during the generation of the rules. So, if the algorithm prim-miner is implemented 
on a flexible platform (e.g., SICStus Prolog with CHR support), when a valid rule 
is extracted it can be immediately simplified with respect to the already generated 
rules. Then, if it is not redundant it can be incorporated at runtime and thus be 
used actively to speed up further resolution and constraint solving steps called by 
the generation algorithm itself. 

Other performance issues. It is likely that the same constraints and their negations 
are candidates for both the Ihs and rhs of the rules. In this case the two following 
optimizations can be used: 

1. If a rule C ^ {d} is generated and not(d) G Candihs then there is no need to 
generate the rule C U {not{d)} false since this rule is trivially redundant. 

2. If a rule Ci =^ {d} is generated and d is also in Candihs then there is no need 
to consider any C2 such that C2 is a superset of Ci containing d to form the 
Ihs of another rule. 
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3. If a rule C U {di} =J> {^2} is produced using the evaluation of the goal 
CU{di, not{d2)} (which fails), and not(d2) £ Candius and not{di) G Candrhs 
then the same goal evaluation is also needed to produce the rule CU{not{d2)} ^ 
{not{di)} . Thus, we can avoid for such constraints this kind of repeated goal 
evaluations. 

For example, the execution of the goal min{X,Y, Z), Y<X, Zj^Y may lead 
to the generation of the following rules: 

min{X,Y,Z), Y<X => Z=Y. 
min{X,Y,Z), Z^Y Y>X. 

2.4 Generation of Primitive Splitting Rules 

Splitting rules have been shown to be interesting in constraint solving ( |Apt 199'8| ), 
since they can be used to detect early alternative labeling cases or alternative 
solution sets. These rules are handled by CHR^ an extension of CHR proposed 
in l|Abdennadher and Schiitz 199811 . Such rules that can be generated by the exten- 
sion proposed in this section are for example: 

and{X,Y,Z), Z^O X=0 V Y=0. 

min{X,Y,Z) X=Z V Y=Z. 

where and{X, Y, Z) means that Z is the Boolean conjunction of X and Y, 
defined by the facts anc?(0, 0, 0), and(l, 0, 0), anc?(0, 1, 0), and(l,l,l); and where 
min{X, Y, Z) is defined by the constraint logic program of Example 12. II 

In the following, we restrict oursclf to the generation of primitive splitting rules. 

Definition 2.4 

A primitive splitting rule is a rule of the form C ^ di\/d2, where di, d2 are primitive 
constraints and C is a set of primitive and user-defined constraints. V is interpreted 
like the standard disjunction, and primitive splitting rules have a straightforward 
associated semantics. 

The modification of the basic prim-MINER algorithm is as follows. For all Cihs, 
it must also consider each different set {(ii,c?2} ^ Crhs with di ^ d2 and check if 
the goal {Baseihs U Cihs U {not{di) , not{d2)}) fails with respect to the constraint 
logic program P. If this is the case, it simply adds to TZ the primitive splitting rule 

(Baseihs U Cihs di V ^2)- 

For example, the rule min{X, Y, Z) ^ X—Z V Y=Z can be obtained by running 
the goal min{X,Y, Z), Xj^Z, Y^Z and checking that its execution finitely fails. 

Here again, the soundness of the generation relies on the properties of the under- 
lying resolution used. 

In practice, a huge number of splitting rules are not interesting because they are 
redundant with respect to some primitive propagation rules. So, to remove these 
uninteresting rules, the generation should preferably be done, by first using prim- 
MINER to obtain only primitive propagation rules, and then using the modified 
version of the algorithm to extract primitive splitting rules. Thus in this second 
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step, redundant rules, with respect to the set resulting from the first step, can be 
discarded on-the-fly. For example the rule a, 5, c => d V e will be removed if we 
have already generated the rule a,b => d. Moreover, in such trivial redundancy 
cases, the test of the validity of the rule a,b,c ^ d V e can itself be avoided, and 
more generally the test of any rule of the form C => d V / and C / V d, where 
{a,b} C C and / G Candrhs- 

3 Generation of More General Propagation Rules 

In this section, we modify the algorithm presented in Section |21 to handle a broader 
class of rules called general propagation rules, that encompasses the primitive prop- 
agation rules, and that can even represent recursive rules over user-defined con- 
straints (i.e., rules where the same user-defined constraint predicate appears in the 
Ihs and rhs). 

Definition 3.1 

A general propagation rule is a failure rule or a rule of the form Ci C2, where 
Ci and C2 are sets of primitive and user-defined constraints. 

The notion of validity defined for primitive propagation rules also applies to this 
kind of rules. 

In the PRIM-MINER algorithm, the validity test of a primitive propagation rule 
Baseihs U Cihs {d} is performed by checking that the goal Baseihs U Cihs U 
{not(d)} fails. For general propagation rules, d is no longer a primitive constraint 
but may be defined by a constraint logic program. In this case, the evaluation should 
be done using a more general resolution procedure to handle negated subgoals. 
However, to avoid the well known problems related to the presence of negation we 
can simply use a different validity test based on the following theorem. 

Theorem 3.1 

Let Ci C2 be a general propagation rule and V be the variables occurring in 
Ci. Let Si be the set of answers {ai, . . . ,a„} to the goal Ci, and 5*2 be the set 
of answers {61, . . . to the goal Ci U C2. Then the rule Ci C2 is valid if 
P\CT h -(3-v((ai V ... V a„) A -3_v(6i V ... V 5„))). 

This straightforward property comes from the completeness result of standard 
CLP schemes l|,TafFar and Maher 1994|l which ensures that if a goal G has a finite 
computation tree, with answers Ci, . . . , c„ then P* , CT \= G ^ ^-Vg (ci V . . . V c„), 
where Vg is the set of variables appearing in G. 

So, the modification proposed in this section consists simply in the replacement 
in algorithm prim-MINER of the call to the goal Baseihs U Cihs U {not{d)} to check 
the validity of the rule Baseihs U Cihs =^ {d}, by the following steps. 

• First, collect the set of answers {ai, . . . , a„} to the goal Baseihs U Cihs- 

• Then, collect the set of answers {61, ... , 6^} to the goal Baseihs U Cihs U {d}. 

• In each answer Oj (resp. rename with a fresh variable any variable that is 
not in Baseihs U Ciha (resp. Baseihs U Cihs U {d}). 
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• Finally, perform a satisfiability test of (ai V . . . V a„) A -i(6i V . . . V 6m). 

• If this test fails then the rule Baseihs U Cihs {d} is valid. 

Since answers only contain primitive constraints and since the set of primitive 
constraints is closed under negation, then we can perform the satisfiability test by 
rewriting (ai V ... V a„) A -1(61 V ... V bm) into an equivalent disjunctive normal 
form, and then use the solver for primitive constraints on each sub-conjunctions. 

It should be noticed that in cases where the evaluation of one of the two goals does 
not terminate before the bound of resolution depth is reached then the propagation 
rule Baseihs U Cihs ^ {d} is not considered as valid and the next rule is processed. 

Example 3.1 

Consider the following user-defined Boolean constraints: neg{X^ Y) imposing that 
Y is the Boolean complement of X and xor{X^ Y, Z) stating that Z is the result of 
the exclusive Boolean disjunction of X and Y. These two constraints are defined by 
a straightforward constraint logic program. Among others, the modified algorithm 
presented above can generate the following rules: 

xor{X,Y,Z), Z=l neg{X,Y). 
xor{X,Y,Z), Y=l => neg{X,Z). 
xor{X,Y,Z), X=l => neg{Y,Z). 

For example, to test the validity of the first rule, the process is the following. First, 
the answers Ai := X=l A Y=0 A Z^l and A2 := X=0 A Y=l A Z=l to the goal 
xor{X, Y, Z), Z=l are computed. Then for the goal xor{X, Y, Z), Z=l, neg{X, Y) 
the same answers are collected. Finally the satisfiability test of {Ai\/ A2) A^{Ai\/ A2) 
fails and thus establishes the validity of the rule. 

One other general rule that can be generated using the modified algorithm pre- 
sented is for example: 

and{X,Y,Z) ^ min{X,Y,Z). 

Furthermore, rules representing symmetries can be automatically detected using 
the modified algorithm. For example, the rules 

min{X,Y,Z) min{Y,X,Z). 
xor{X,Y,Z) ^ xor{Y,X,Z). 

expressing the symmetry of the 771m and the xor constraints with respect to the 
first and second arguments can be generated. In ( |Abdennadher and Rigotti 2001a| ), 
it has been shown that these rules are very useful to reduce the size of a set of prop- 
agation rules since many rules become redundant when we know such symmetries. 

However, it must be pointed out that if the test presented in this section allows 
us to handle a syntactically wider class of rules, it relies on different goal calls than 
the test of Section ITTI So when testing the validity of a primitive propagation rule 
one of the techniques may lead to terminating evaluation while the other one may 
not. Thus in the case of primitive propagation rule it is preferable not to replace 
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one test by the other, but to use both m a complementary way (run one of them, 
and if it reaches the bound of resolution depth then apply the other) . 

4 Generation of Simplification Rules 

Since a propagation rule does not rewrite constraints but adds new ones, the con- 
straint store may contain superfluous information. Constraints can be removed from 
the constraint store using simplification rules. 

Definition 4-1 

A simplification rule is a rule of the form C'l <^4> C2, where Ci and C2 are sets of 
primitive and user-defined constraints. 

In this section, we show how the rule validity test used in Section |31 can be 
applied to transform some propagation rules into simplification rules. For a valid 
propagation rule of the form C ^ D, wc try to find a proper subset E oi C such 
that DU E ^ C is valid too. If such E can be found, the propagation rule C ^ D 
can be transformed into a simplification rule of the form C D U E. 

To simplify the presentation, we present an algorithm to transform (when pos- 
sible) propagation rules into simplification rules independently from the algorithm 
presented in Section |3| Note that the algorithm for the generation of propagation 
rules can be slightly modified to incorporate this step and to directly generate 
simplification rules. 

The algorithm is given in Figure [5] and takes as input the set of generated prop- 
agation rules and the common part that must appear in the Ihs of all rules, i.e. 
Baseihs- 



begin 

P' ■- P 

for each propagation rule of the form C => D in P do 

Find a proper subset E oi C such that Baseihs 2 ^ ^-nd 
D U -E => C is valid (using the validity test of Section I^J 
If E exists then 

P' ■- {P'\{C ^ D})VJ{C ^ D\JE} 
endif 

endfor 

output P' 

end 



Fig. 2. The Transformation Algorithm 

To achieve a form of minimality based on the number of constraints, we generate 
simplification rules that will remove the greatest number of constraints. So, when we 
try to transform a propagation rule into a simplification rule of the form C <^ DUE 
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we choose the smaUest set E (with respect to the number of atomic constraints in 
E) for which the condition holds. If such a i? is not unique, we choose any one 
among the smahest. The conchtion Baseihs % E needed to be able to transform 
the propagation rules into simplification rules that rewrite constraints to simpler 
ones (primitive constraints if possible), as shown in the following example. 

Example 4-1 

Consider the following propagation rule R generated for the append constraint with 
Baseihs = {append{X, Y, Z)}: 

appendix, Y,Z), X^[] Y^Z. 

The algorithm tries the following transformations: 

1. First, it checks if rule R can be transformed into the simplification rule 
append{X,Y, Z), X=[] <^ Y=Z. This is done by testing whether the rule 
Y=Z append{X,Y, Z), X=[] is valid. But, this is not the case and thus 
the transformation is not possible. 

2. Next, the algorithm finds out that the rule Y—Z, X^[] => append{X,Y, Z) 
is a valid rule and then the propagation rule R is transformed into the sim- 
plification rule append(X, F, Z), X=[]^X=[], Y=Z. 

Note that the propagation rule Y=Z, append{X, Y, Z) ^ -'^'=0 is also valid but 
transforming rule R into 

appendix, Y,Z), ^ append{X,Y, Z), Y=Z. 

will lead to a simplification rule which is uninteresting for constraint solving, i.e. 
using this rule the append constraint cannot be simplified and remains in the con- 
straint store. The algorithm disables such transformation by checking the condition 
that all constraints of Baseihs are not shifted to the right hand side of the rule 
{Baseihs % E). 

5 Implementation of the Generated Rules in CHR 

The generated rules may contain constraints that are built-in constraints for the 
CHR system. To have a running CHR solver, these constraints have to be en- 
coded in a specific way. First, equality constraints appearing in the left hand side 
of a rule are propagated all over the constraints in its left and right hand side. 
Then the resulting constraints are simplified. This can be performed as follows. 
In turn each equality constraint appearing in the Ihs is removed and transformed 
in a substitution that is applied to the Ihs and the rhs. Then the next equality 
constraint is processed. For example, the simplification rule and{X,Y, Z), Z=l <^ 
X=l, Y=l, Z=l will be transformed into and{X,Y,l) <^ X^l, Y=l. Secondly, 
for other built-in constraints the transformation leads to CHR rules containing a 
guard l|Friihwirth 1998|l . For example, if < is a built-in constraint of the CHR sys- 
tem then the rule min{X,Y, Z), X<Y <^ Z=X, X<Y is transformed into the 
guarded CHR rule mm(X, F, Z) ^ X<Y \ Z^X. 



14 



S. Abdennadher and C. Rigotti 



6 Related Work 

In | |Abdennadher and Rigotti 2001a| ), a method has been proposed to gener- 
ate propagation rules from the intentional definition of the constraint predicates 
(eventually over infinite domains) given by mean of a constraint logic program. 
It extended previous work ( |Apt and Monfroy 1999| |Ringeissen and Monfroy 2000| 
[Abdennadher and Rigotti 2000| ) where different methods dedicated to the genera- 
tion of propagation rules for constraints defined extensionally over finite domains 
have been proposed. Compared to the work described in (l Abdennadher and Rigotti 2001a| ) 
the approach presented in this paper has several advantages: 

• It enables user-defined constraints to occur in the right hand side of rules, 
while this is not handled by ( [Abdennadher and Rigotti 2001a| ). As a by- 
product rules representing symmetries as 

min{X, Y, Z) min{Y, X, Z). 

can then be automatically detected. 

• It allows the generation of splitting rules like 

append{X,Y,Z), Z=[A] X=[A] V Y=[A] 

supported by the extension of CHR called CHR^ ijAbdennadher a nd Schiit z 1998|l . 
These rules have been shown to be interesting in constraint solving ( |Apt 1998| ), 
but have not been considered in ( |Apt and Monfroy 1999||Ringeissen and Monfroy 2000| 
[Abdennadher and Rigotti 2000| [Abdennadher and Rigotti 2001a| ). 

• Even if we restrict our attention to the class of rules handled in ( [Abdennadher and Rigotti 2001al ) 
the approach presented in this paper leads to a more expressive set of rules. 

For example, the rule 

appendix, Y,Z), Y=[] X=Z. 

cannot be generated by the approach proposed in ( [Abdennadher and Rigotti 2001a[ | 
while the algorithm described in Section |21 is able to obtain it by executing 
the goal append{X, Y, Z), Y=[], Xy^Z with a tabled resolution for CLP. 

• Additionally, it needs less information about the semantics of the primitive 
constraints. For example it generates the rule 

min{X,Y,Z) =^ Z<X, Z<Y. 

simply by calling the goals min{X,Y, Z), Z>X and min{X,Y, Z), Z>Y, 
and checking that their executions fail. While in this case, the algorithm pre- 
sented in ( [Abdennadher and Rigotti 2001aj | requires some extra information 
concerning the semantics of < (information that would be provided in general 
by the user). 

In this paper, we also described a method to transform propagation rules into 
simplification rules. It has been shown in ( [Abdennadher and Rigotti 2001b[ ) that 
this transformation has a very important impact on the efficiency of the solver pro- 
duced. In ( [Abdennadher and Rigotti 2001b[ ) the propagation rules are modified to 
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obtain (when possible) simplification rules using a technique based on a confluence 
notion. This is a syntactical criterion that works when we have at hand the whole 
set of rules defining the constraint. Thus it cannot be applied safely if only a part 
of the propagation rules have been generated. It also requires a termination test 
for rule-based programs consisting of propagation and simplification rules, and this 
test is undecidable for some classes of programs. The new transformation method 
presented in this paper avoids these two restrictions. 

The generation of propagation and simplification rules is also related in some as- 
pects to Generalized Constraint Propagation IjLe Pro vost and Wallace 1993(1. Gon- 
structive Disjunction ( |Van Hen tenryc k et al. 1998| iWiArtz and Miiller 1996| . and 
Inductive Logic Programming ( Mugglct on and De Raedt 1994| ) as briefly discussed 
in ) |Abdennadher and Rigotti~2001a| ). However, it should be pointed out that to our 
knowledge these works have not been used for the generation of constraint solvers. 

7 Conclusion and Future Work 

In this paper, we have presented three algorithms that can be integrated to build 
an environment to help solver developers when writing CHR programs. 

The approach described allows the generation of CHR propagation and simplifi- 
cation rules. It can be applied on constraints defined over finite and infinite domains 
by mean of a constraint logic program. Moreover, it enables the developer to search 
for rules having such user-defined constraints in both their left and right hand sides. 

We have also shown that compared to the algorithms described in ( |Apt and Monfroy 1999| 
IRingeissen and Monfroy 2000||Abdennadher and Rigotti 2000||Abdennadher and Rigotti 2001al ) 
to generate rule-based constraint solvers, this approach is able to generate more ex- 
pressive rules (including recursive and splitting rules). 

One interesting direction for future work is to investigate the integration of con- 
structive negation (e.g., ( |Stuckey 1995| )) in tabled resolution for CLP to generate 
constraint solvers, in order to check the validity of the propagation and simplifica- 
tion rules in more general cases. Another complementary aspect is the completeness 
of the solvers generated. It is clear that in general this property cannot be guaran- 
teed, but in some cases it should be possible to check it, or at least to characterize 
the kind of consistency the solver can ensure. 
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